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^ ' Abstract. The flow characterized by a linear longitudinal velocity field Ux{x,t) = a{t)x, where a(f) = ao/(l a uniform 

, density n{t) o= a(?), and a uniform temperature T{t) is analyzed for dilute granular gases by means of a BGK-like model 

kinetic equation in d dimensions. For a given value of the coefficient of normal restitution a, the relevant control parameter 
^ ' of the problem is the reduced deformation rate a*{t) = a{t)/v{t) (which plays the role of the Knudsen number), where 

, v(f) °^ n{t) \/T{t) is an effective collision frequency. The relevant response parameter is a nonlinear viscosity function r]*(a*) 

defined from the difference between the normal stress Pxx{t) and the hydrostatic pressure p{t) = n{t)T{t). The main results 
of the paper are: (a) an exact first-order ordinary differential equation for T]*(a*) is derived from the kinetic model; (b) 
a recursion relation for the coefficients of the Chapman-Enskog expansion of r]*{a*) in powers of a* is obtained; (c) the 
Chapman-Enskog expansion is shown to diverge for elastic collisions (a = 1) and converge for inelastic collisions (a < 1), 
fir-^ . in the latter case with a radius of convergence that increases with inelasticity; (d) a simple approximate analytical solution for 

^ ' 77*(a*), hardly distinguishable from the numerical solution of the differential equation, is constructed. 
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g i A granular gas is a large collection of mesoscopic or macroscopic particles which collide inelastically and are 

Q ' maintained in a fluidized state. The prototypical model of a granular gas consists of a dilute system of smooth inelastic 

hard spheres characterized by a constant coefficient of normal restitution a < 1 [ 1]. A kinetic theory description based 
on the Boltzmann equation for inelastic collisions has proven its relevance for the understanding of the main properties 
of dilute granular gases, such as non-equipartition of energy and high-velocity tails in uniform states or non-Newtonian 
transport properties in non-uniform steady states. 
\ The kinetic theory definition of nonequilibrium temperature T as a measure of the mean kinetic energy (in the 

■ Lagrangian frame) per particle can be directly extended to granular gases (except that it is customary to measure T in 
energy units, i.e., to take the Boltzmann constant equal to unity). The key difference with respect to normal gases is 

Q>^ i that the inelastic character of collisions gives rise to an energy sink, so that the temperature change per unit time due 

' to collisions is <3r7"|coii = ^^T, where C, is the so-called cooling rate. It is given (in the Maxwellian approximation) by 

OO ! ^ 

C = ^(l-a)c., <o^j^^^^^na (1) 

k>( \ where d is the dimensionality of the system, n is the number density, a is the diameter of a sphere, and m is its 

5_j ■ mass. Application of the Chapman-Enskog (CE) method to the inelastic Boltzmann equation allows one to derive the 

I Navier-Stokes (NS) constitutive equations for inelastic hard spheres [2]. In particular, Newton's law reads 

/ dui duj 2„ - \ 
P. = p5,-r,.s(^ + ^~^V.u5,j, (2) 

where f,,- is the pressure tensor, p = {\/d)TxP = nT is the hydrostatic pressure, u is the flow velocity, and tjns is the 
NS shear viscosity. The latter quantity is given, in the simplest Sonine approximation, by |i2|J 



nT 1 + a 



®. (3) 



The uniform longitudinal viscous flow is an unsteady compressible flow defined by a linear longitudinal velocity 
field, a uniform density, and a uniform temperature Istl- The exact balance equations for mass and energy yield 

Ux{x^t) = a{t)x, n{t)~—a{t), a{t) — — , (4) 

ao 1 + dQt 



a>0 



a<0 
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FIGURE 1. Sketch of the uniform longitudinal flow for ao > (left) and ao < (right). 



d'T{t) = -^^P.At)-at)T{t). (5) 

This simple flow is also known as homo-energetic extension and, along with the uniform shear flow P, 's'], is a 
particular case of a more general class of homo-energetic affine flows characterized by d^Ui/dxjdxj^ = |4]. The 
initial longitudinal deformation rate rate oq is the only control parameter determining the departure of the fluid from 
the homogeneous state, thus playing a role similar to that of the shear rate in the uniform shear flow state. On the other 
hand, in contrast to the uniform shear flow, the sign of oq plays a relevant role and defines two distinct situations (see 
Fig. [1]). The case oq > corresponds to a progressively slower expansion of the gas from the plane x = into all of 
space. On the other hand, the case oq <Q corresponds to a progressively faster compression of the gas towards the 
plane x~Q. The latter takes place over afinite time period t = |flo| However, since the collision frequency rapidly 
increases with time, the finite period t = |flo|^' comprises an infinite number of collisions per particle (see below). 
According to Newton's law (|2]i, the normal stress P^x in the uniform longitudinal flow is given by 

l^it) = -— = 1 —2 — — , a n = — -— <:x: , (o) 

p{t) d i + c*/2 ^' v(f) «oyr(0 

where ^* = ^/v is a dimensionless constant (except that it depends on d and a). Note that T* is the ratio between 
the "longitudinal" temperature = Pxx/n and the true granular temperature T . The absolute value of the reduced 
longitudinal rate a* plays the role of the Knudsen number of the problem. The linear relationship (|6]l between 7^ and 
a* is expected to hold for small values of |fl*| only. More in general, one can define a (reduced) nonlinear viscosity 
function r\*{a*) by the relation 

r;(0 = l-2^T7*(fl*(OK(f), Hni^T7*(fl*) = ^-j-i^^co. (7) 

The expansion of 77*(a*) in powers of a* yields the CE series 

^*{a*) = ^Y^c,a*\ (8) 

k=Q 

where c\ is a Burnett coefficient, C2 is a super-Burnett coefficient, and so on. 

In this paper a kinetic theory description of the longitudinal viscous flow defined by Eqs. (01 and (|5]) is presented. 
As will be shown, a steady-state value of the temperature is reached only if aq < 0. The main result is the derivation of 
a nonlinear first-order differential equation for r\*{a*). This equation allows one to get the CE coefficients Cj^ of Eq. (O 
in a recursive way. The results show that the CE series is convergent for inelastic hard spheres, while it diverges in the 
elastic case. The viscosity function Tj * (a* ) for any value of the inelasticity parameter a can be obtained by numerically 
solving the differential equation with appropriate boundary conditions. Finally, a simphfied analytical solution hardly 
distinguishable from the numerical solution is proposed. 



KINETIC MODEL 

In order to obtain explicit results, it is convenient to consider the following BGK-hke kinetic model |0] of the inelastic 
Boltzmann equation: 

(5, + V • V)/ = - V(/ - /he,) + [(v - u)/] , (9) 



where / is the velocity distribution function and /hcs is the local version of the homogeneous cooling state distribution 
111]. The NS viscosity coincides with that of the Boltzmann equation if v is given by Eq. However, a simpler 
and more consistent choice is v = 5(1 + a)co ff], so that ^* = [{d + 2)/2d]{l — a). In the uniform longitudinal flow 
the velocity distribution function /(x, v,f ) becomes spatially uniform when the velocities are referred to a Lagrangian 
frame moving with the flow, i.e., f{x, v, f ) = /(V, t), where V = v — u(jic, f ) is the peculiar velocity. After simple algebra, 
Eq. (|9]l can be rewritten as 

{d, - flo5y,V07 = -V (7- 7hcs) + 1^ • (V7) , (10) 

where 

7(V,T)^^/(x,v,f), v(T)^^v(r), CW^^CW, T^ao'ln(l+«oO- dD 
n[t) n[t) n[t) 

Equation ( fTOl i shows that the original uniform longitudinal flow problem can be mapped onto the equivalent problem 
of a uniform gas with a velocity distribution / and subject to the action of a non-conservative force —maoVxX. The 
density associated with / is the initial density no and thus remains constant. In fact, / is the probability distribution 
function normalized to iiq. The time variable T = /q dt'n{t') /no is a nonlinear measure of time scaled with the number 
density. It is unbounded even if oq < since in that case t — > 00 when f ^ | ao | ' ■ 
Taking velocity moments on both sides of Eq. ( fTOl ) one gets 

d,T + -aQT, = -XT, d,T, + 2aaT,^-v{T,-T)-XT,. (12) 
a 

If flo > 0, dT:T < even in the elastic case (^ = 0). As a consequence, the gas cools down and the reduced longitudinal 
rate a* ~ a/v ^ increases without bound, i.e., limT^oofl*(T) = 00. On the other hand, in the compression 

case (flo < 0) the "viscous heating" term {2/d)\ao\Tx competes with the "inelastic cooling" term ^T, so that the 
temperature either decreases or increases with time depending on whether in the initial state one has > {2/d)\ao\Tx 
or < {2/d) \ao\Tx, respectively. In either case, a steady-state temperature is eventually reached when both effects 
cancel each other Therefore, if aq < one has limT^oo a* (t) = a* . The steady-state values for the reduced longitudinal 
rate a*, the temperature ratio T^/T, and the reduced viscosity 77* can be easily obtained from Eq. (fT2] i and are given 
by 

a = C — T = — n = — (13) 

2^ 1+dC*' 1 + C*' (1 + C*P' 

Comparison between 77* and the NS value cq in Eq. Q shows that the steady state is inherently non-Newtonian, even 
in the quasi-elastic limit ^* ^ 1, in which case co w 1 — C*/2 whilst Tj * w 1 + {d — 2)^*. Note that in the elastic case 
(1^ = 0) with flo < only the viscous heating term is present in the energy balance equation, so that limT^oc a* (t) = 0. 

In summary, the zero Knudsen number value (a* = 0) is a "repeller" of the time evolution of fl*(T) in the inelastic 
case. For elastic collisions, the state a* = is also a repeller if uq > but it is an "attractor" of a*(T) if oq < 0. Figure 
|2|a) shows a*, T*^, and rj* as functions of a for three-dimensional systems (d = 3). It is interesting to note that the 
steady-state reduced viscosity 77* is nearly independent of a. 

As noted before, if oq < the whole compression process takes place during a finite time period f = | aq | ^ ' . Let us see 
now that the corresponding number of collisions per particle is infinite. The number of collisions per particle during an 
elementary time interval df is ds — vdt = vdz. The accumulated number of collisions between the initial state and the 
scaled time T is then ^(t) ~ Vq Jq dz' ^T{t')/Tq, where Vo and Tq are the initial collision frequency and temperature, 
respectively. If oo < the temperature is bounded between Tq and Ts, i.e., min(ro, r, ) <T{t') < max(ro, 7i). Therefore, 
min(l, ^JTs/Tq)Vq'Z < s{z) < max(l, ^/^]/^}))Vo^:■ Since t ^ 0° when t |flo| ^ it follows that 5 ^ 00 in that limit. 



NONLINEAR VISCOSITY 

In order to get the whole rheological function 77*(a*) we need to eliminate time in favor of a*(T) = flo/v(T), taking 
into account that both v and ^ are proportional to [^(t)]'', where q = 5. It is convenient to consider for the moment q 
as an arbitrary positive parameter and setq=j at the end. Thus, since — T* T, one has 



T:Tjc = T^dtT + Ti^d-ca )da' = [T^ - qa da* ) d^T. 



(14) 
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FIGURE 2. (a) Plot of the steady-sate values a*, T*^, and T]* as functions of a. (b) Ratio |c^/c^+i| for a = 0.5, 0.9, and 1; the 
horizontal dashed lines indicate the corresponding values of (c) Nonlinear viscosity for a = 0.5 and 1 as obtained from the 
numerical solution of Eq. J15b (solid lines) and from Eq. i20\ (dash-dotted lines); the circles denote the steady-state points and the 
arrows indicate the direction of the time evolution. All the panels correspond to three-dimensional hard spheres (q = ^, d = 3). 



Making use of Eq. (fT2] l it is easy to get the following nonlinear first-order differential equation for the function T7*(fl*): 

iri*+a*da,ri*)=0. (15) 
Inserting the CE expansion ([8]) one can get the coefficients q in a recursive way: 



a \ a 



Figure |2|b) shows the absolute value of the ratio c^^/c^+i for three-dimensional hard spheres (q = ^, d — 3) with 
a — 0.5 (highly inelastic case), a — 0.9 (moderately inelastic case), and a — I (elastic case). In the inelastic cases 
the ratio Iq/c^^+i | tends to a non-zero finite value, thus indicating that the CE series ^ converges. In fact, the results 
are consistent with lim^-^cx, | = i.e., the radius of convergence of the CE series coincides with the steady- 

state value of the Knudsen number. In contrast, in the elastic case one has lim<.^oo Iq/cj^+i | = 0, so that the radius of 
convergence shrinks to zero and the CE series ( fT6b diverges. In fact, if ^* = 0, the magnitude of the coefficients cj^ 
grows so fast that the summation term in Eq. ( fTSl ) can be neglected and one has Iq/c^-^i | « {d/2q)k^^ . 
Inspection of Eq. ( fTSl) shows the following asymptotic behaviors in the limit of jo* | oo-. 



flo > : lim«*^„oTj*(fl*)fl* ^ lim„.^oor;(fl*) = 0, 

flo<0: lima*^_«,T7*(fl*)fl* lim„*^_oo r/(fl*) = 



(17) 



Therefore, in the limit a* ^ oa there is no motion (relative to the flow velocity) along the x direction, while in the 
opposite limit a* — oo the motion becomes one-dimensional along the x axis. 

The complete viscosity function 77*(fl*) for any a* , including values of the Knudsen number beyond the radius 
of convergence |a*|, is obtained by numerically solving the differential equation ( fTSl l with appropriate boundary 
conditions. To get the branches a* > and a* < a* < one must start from a small value of |fl*| and apply the 
initial condition r\* = cq + c\a* . The branch a* < a* < is obtained by starting from a large value of a* and applying 
the boundary condition Tj* — —d/2a*. The resulting curve r\*{a*) is plotted in Fig.|2c) for three-dimensional hard 
spheres (^ = j, c/ = 3) with a = 0.5 and 1. The arrows on the curves indicate the direction followed by the time 
evolution of a*{x). 

While the differential equation ( fTSl l can be solved numerically, it is interesting to approximate its solution by an 
analytical expression. To that end, let us see ^ as a small perturbation parameter and formally expand the solution to 
Eq. ( fTSb in powers of q: 

TJ*(fl*) = 77o*(fl*)[l-^/!l(fl*) + ---]. (18) 



FIGURE 3. Surface plot showing the dependence on both a* and a of of the viscosity function rj* given by Eq. i20\ . The thin 
and thick lines refer to the NS and steady-state points, respectively. 



By inserting Eq. ( flSl l into Eq. ( fTSl ) and keeping terms linear in q, one gets, after some algebra. 



T7o*(«*) = 



d-2 



Pic 



hi{a*) = 



Pic 



4j32(fl*) 




Note that the zeroth-order function T7o(fl*) is independent of the reduced cooling rate ^* 
higher order terms in the ^-expansion ( fTSl ). let us consider the simple Pade approximant 



1 +qh\{a* 



(19) 

Instead of proceeding to 
(20) 



This approximation has the merit of providing correctly (for any value of q) the values of the NS coefficient cq [see 
Eq. ( fT6b l and the steady-state point r\* — r\*{a*) [see Eq. (fT3])1. as well as the limiting behaviors ( fTTb . Moreover, even 
in the physical case of hard spheres {q = 5), Eq. ( |20l ) represents an excellent analytical approximation to the numerical 
solution of Eq. ( fTSl ). especially for highly inelastic systems, as shown in Fig.|2fc). The discrepancies are only apparent 
near the maximum of 77* (a*) and in the decay to the left of the maximum. Figure |3] shows a surface plot of Eq. ( l20b . 



CONCLUSIONS 

The uniform longitudinal viscous flow is an unsteady compressible flow [see Eqs. (|4|i and Q] that, despite its apparent 
simplicity, constitutes a non-trivial playground for nonequilibrium statistical mechanics beyond the NS description. 
While this state has been studied in the past for normal gases 1 3], the present work addresses the problem in the case of a 
granular gas of inelastic hard spheres. The relevant control parameter is the reduced longitudinal rate a*{t)= a{t)/v{t) 
(whose magnitude plays the role of the Knudsen number of the problem) and the relevant response function is the time- 
dependent generalized viscosity 77 * (a* ) defined by Eq. (|7]i. The exact energy balance equation ^ shows that if aq > 
(expansion states), the temperature monotonically decreases and so a* (t) increases with time, for both elastic (^ = 0) 
and inelastic (^ > 0) collisions. On the other hand, if aq < (compression states), a* (f ) goes to zero from below in the 
elastic case, whereas it tends to a stationary value a* < in the inelastic case [see Fig.|2a)]. Expressed in other terms, 
one can say that the homogeneous state (aq = or a* = 0) is unstable against any perturbation oq 0, no matter how 
weak it is, if the particles are inelastic. On the other hand, a* = is stable for elastic collisions if ao < 0. Given that 
the CE expansion dHJ is carried out about (and measures the departure from) the reference homogeneous state (a* = 0), 
it follows the "arrow of time" for granular gases but goes against it for normal gases (if ao < 0). As a consequence, it 
can be expected on physical grounds that the CE expansion diverges in the case of elastic colhsions but does converge 
if the collisions are inelastic, the radius of convergence |fl*| increasing with inelasticity. 



To confirm the above expectation and also to get the rheological function Tj* (a*), the simple model kinetic equation 
^ has been considered in this paper. The viscosity function J]* (a*) predicted by the model is the solution to the 
nonlinear first-order differential equation (fTsT i (with q = 5). The CE coefficients c^. are shown to verify the asymptotic 
law limit^oo Iq/q+iI = \a*\, thus confirming the convergence of the series ([8]l with a finite radius a* — \a*\, which 
vanishes in the elastic limit [see Fig.|3b)]. The function rj*{a*) exhibits an interesting behavior [see Fig.|2tc)]. Starting 
from the NS value 77* (0) = cq, it decreases as the longitudinal rate a* increases, behaving as 77* (a*) « [d /2{d— !)]/«* 
in the limit of large positive a* . In the latter limit, the gas becomes highly anisotropic with no random motion along the 
flow direction. The situation is even more interesting in the case of negative a* . As the magnitude of the longitudinal 
rate increases, the viscosity function starts increasing, reaches a maximum value past the steady-state point, and 
decreases thereafter, behaving as Tj*(fl*) sa {d/2)/\a*\ in the limit of large The gas becomes highly anisotropic 
also in this limit, but this time because the degrees of freedom orthogonal to the flow direction are suppressed. As Fig. 
|2|c) shows, the largest quantitative influence of inelasticity on the shape of rj*{a*) occurs near the maximum. It is also 
in this region where the analytical approximation (|20] | exhibits the largest deviations from the numerical solution of 
Eq. (fTsT l. Otherwise, Eq. ( |20] | succeeds in capturing the main features of the solution, including the steady-state and 
NS values, as well as the asymptotic behaviors. However, an important property not accounted for by Eq. JTSl l is the 
divergence of the CE series in the elastic case (C* = 0). This is a consequence of the fact that Eq. ( fTST i is based on the 
^'-expansion ( fTST i (truncated after the first term) and the CE expansion in the reference system (q ~ 0) is convergent, 
regardless the value of ^* . 

As a final point, it is worthwhile rewriting the CE expansion (O in the dimensional form 

d—1 ^ Ci. / dur\^^^ 
P« = p-2-— 77nsE-i^K^ . (21) 

When arbitrary hydrodynamic gradients are present, Eq. ( |2TI ) represents a partial contribution to the general CE 
expansion of the normal stress P^^ . The advantage of the uniform longitudinal flow is that it possesses an only 
hydrodynamic gradient, du^/dx, which is spatially uniform. Therefore, in this case the general CE expansion of Pxx 
reduces to Eq. (I2TI 1. However, the convergence or divergence of the partial series (l2Tl i does not depend on whether the 
system is actually in the uniform longitudinal flow or in any other state. According to the results derived in this paper, 
it turns out that the partial series (1211 1 diverges for normal gases but converges for granular gases. A similar conclusion 
is reached in the case of the expansion of the shear stress Pxy in powers of the shear rate dux/dy 
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